Heating Effects in a Chain of Quantum Dots 
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We study heating effects in a chain of weakly coupled grains due to electron-hole pair creation. 
The main mechanism for the latter at low temperatures is due to inelastic electron cotunneling 
processes in the array. We develop a quantitative kinetic theory for these systems and calculate 
the array temperature profile as a function of grain parameters, bias voltage or current, and time 
and show that for nanoscale size grains the heating effects are pronounced and easily measurable 
in experiments. In the low- and high-voltage limits we solve the stationary heat-flux equation 
analytically. We demonstrate the over-heating hysteresis in the large-current or voltage regimes. 
In addition we consider the influence of a substrate on the system which acts as a heat sink. We 
show that nano dot chains can be used as highly sensitive thermometers over a broad range of 
temperatures. 

PACS numbers: 72.15.Jf, 73.63.-b, 85.80.Fi 



I. INTRODUCTION 

Great efforts in contemporary materials scienee re- 
search focus on the understanding of thermal properties 
of nanoscale deviccsEfQ. The electron transport through 
arrays of nanogranular systems is well understoodErt3 at 
small currents when the charge carriers have enough time 
to equilibrate their temperature with the temperature of 
the phonon bath. It is known that most of the trans- 
port characteristics of the system, e.g., the resistance, 
are in fact very sensitive to the temperature of the elec- 
tronic excitations in the granules. It is much less known 
about self-heating effects in these systems when the cur- 
rent is not very small and the charge flow makes the 
internal temperature different (higher) than the temper- 
ature of the phonon bath. Understanding of the heating 
effects on properties of nanodevices is especially impor- 
tant for practicaLaonlications. In addition, recent exper- 
imental researchlixtS has focused on the application of 
the thermal properties of low-dimensional devices for ef- 
ficient power generation. Granular thermoelectric mate- 
rials have the advantage that one can control the system 
parameters and therefore the device properties. Due to 
the high sensitivity of the resistance of granular materials 
on temperature these materials are good candidates for 
broad band temperature sensors that can operate with 
high resolution from cryogenic up to room temperatures. 
This defines an urgent quest for a quantitative descrip- 
tion of heating effects on properties of nanodevices. 

In this paper we investigate the effect of inelastic co- 
tunneling on heating of weakly coupled nano dot chain 
of length N coupled to two leads (see Fig. nl). We calcu- 
late the grain temperature profile as a function of grain 
parameters, bias current, and time for different system 
setups of one-dimensional chains, which may or may not 
be coupled to a substrate through electron-phonon colli- 



Each grain is characterized by two energy scales: (i) 
the mean energy level spacing 6 and (ii) the charging 
energy Ec = e^(47rKeoa)~^, where k is the relative static 
permittivity of the grain material. We concentrate on the 
case of metallic grains which are most commonly encoun- 
tered in applications satisfying the inequality Ec ^ 5. 

In addition to the above energy scales the system under 
consideration is characterized by the tunneling conduc- 
tances pj . In this paper we concentrate on the regime of 
Coulomb blockade, meaning that the tunneling conduc- 
tances are much smaller than the quantum conductance, 
5( ^ 1. Our considerations are valid for temperatures 
6<T<Ec. 
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FIG. 1. (Color online) Sketch of a weakly coupled nano dot 
chain of length A'" connected to two leads, with tunneling con- 
ductances (j( <C 1 (where i ± 1/2 denotes the tunneling 
barriers right/left of grain i) . The leads have the grain indices 
and N + 1. The lead temperatures Ti and T2 and substrate 
temperature Ts are fixed. The grain temperature Tg de- 
pends on these temperatures and on tlie applied current or 
voltage. The main mechanism of heating in the grain due 
to inelastic cotunneling processes is depicted by the electron- 
hole pairs inside the grain which have a typical energy of Tg . 
The colors of the grains show a typical distribution of grain 
temperatures throughout the system: the coldest grains are 
close to the leads and the hottest in the middle of the chain. 



II. MECHANISM FOR GRAIN HEATING 

In the considered low temperature region, the only 
transport mechanism contributing to grain heating is 
due to inelastic cotunneling. Cotunneling allows simul- 
taneous charge transport through several junctions by 
means of cooperative electron motion. Single cotunnel- 
ing ffiig- 0), introduced in Ref. Q, see also the re- 
viewE£l. provides a conduction channel at low applied bi- 
ases, where otherwise the Coulomb blockade arising from 
electron-electron repulsion would suppress the current 
flow. The essence of a cotunneling process is that an 
electron tunnels via virtual states in intermediate gran- 
ules thus bypassing the huge Coulomb barrier. This can 
be visualized as coherent superposition of two events: 
tunneling of the electron into a granule and the simul- 
taneous escape of another electron from the same gran- 
ule. There are two distinct mechanisms of cotunneling 
processes, elastic- and inelastic cotunneling. Elastic co- 
tunneling means that the electron that leaves the dot has 
the same energy as the incoming one. In the event of in- 
elastic cotunneling, the electron coming out of the dot 
has a different energy than the entering electron. This 
energy difference is absorbed by an electron-hole excita- 
tion in the dot, which is left behind in the course of the 
inelastic cotunneling process, depicted in Fig. |l|. Below 
we concentrate on inelastic cotunneling because only this 
transport mechanism contributes to heating effects. In 
particular, elastic cotunneling and sequential tunneling 
do not create electron-hole pairs in the grain. 

For the quantitative description of the non-equilibrium 
effects leading to the heating of the chain, we consider 
(i) the heating due to electric resistance, (ii) the heat 
removal due to electron and phonon thermal conductivi- 
ties, and (iii) in some cases the electron-phonon coupling 
of the array to a substrate. The electric conductivity and 
the electronic part of the thermal conductivity are gov- 
erned by inelastic cotunneling processes. Using all these 
ingredients we can write down the kinetic equation for 
the grain temperatures in the following section, describ- 
ing the time evolution and final temperature profile over 
the chain. 



III. KINETIC HEAT EQUATION 

For the temperatures of each grain i in the array we 
can write down the following kinetic equation 
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where ch = (l/3)fc^Tg V is the heat capacitance of each 
grain with ly = l/{a^S) being the density of states in the 

grain; dTg /dt describes the change of the grain tem- 
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peratures Tg in time t. a\^ is the cotunneling conduc- 
tance (defined below) for grain «, V'^^^ being the voltage 



drop and a the grain size; k'^^^ is the thermal conduc 
tivity, with k*-'-* 
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and kV the phonon contribution to the thermal conduc- 
tivity: for a single metallic grain weakly coupled to its 
neighbor ji^^'pt^^ Jcads these given by the following ex- 
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In Eq. (^) and (|2q) 7e/ph are numerical coefficients 
(7e = 327r3/15 and 7ph = Stt^/IS, see 0), /ph the 
phonon mean free path Zph ~ a, and Qd the Debye tem- 
perature. Notice, that the temperature entering Eqs. (0) 
is the grain temperature, Tg , which determines the 
number of typically involved internal energy states of the 
grain. The tunneling conductances gi are to be under- 
stood as the averaged quantities y gt St i with 
i ± 1/2 denoting the tunneling barriers right/left of grain 
i. 

The symbol V in Eq. (nl) refers to the discrete gradient 
along the chain. The last term in the r.h.s. of Eq. (|^), 
Vs describes the heat dissipation from grain i due to 
electron-phonon coupling of the array to a substrate and 
is given byH 
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with Tg being the substrate temperature and E^ the 
electron-phonon coupling constant, which is of the or- 
der 10^W/(m'^K'''). This electron-phonon term is written 
for the clean case meaning that the phonon mean free 
path Xpfi is smaller than the electron length, Xph < lei,' ■ 
Here, we notice that the applied current density j does 
not influence the substrate term. 

Equation (1^) has a transparent physical meaning: The 
energy of the electron-hole pair is released as heat, see 
also Fig. |l| whereas the thermal conduction to the leads 
removes the heat from the grain array until a steady state 
is reached. The typical energy of electron-hole pairs is of 
order of the grain temperatures, Tg , and included in the 

expressions for a^^ and k[* . The leads are assumed to 
be a heat sink and stay at ambient temperature T since 
they are much larger than the grains. The first term on 
the r. h. s. of Eq. (|l|) describes the energy released in the 
grain from an electron-hole pair times the square of the 
discrete gradient of voltage {V^'^'' /a). The latter ensures 
that the result is invariant under a sign change of the 

voltage y(*). The second term V (Ac'^'wTg J describes 



the heat- flux out of grain i. Our eonsiderations are vaUd 
as long as the applied voltage does not break down the 
Coulomb blockade, i.e. when eT^^'^ < Ec- 

In typical experiments the current is kept constant and 
Eq. (|l|) can be written as 
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B. Discrete kinetic heat equation 

Here we present the explicit expression of the kinetic 
heat Eq. (H) for the discrete chain of grains. First, 
we derive an expression for the dimensionless voltage 
drop on grain (j), which follows from V^*-* = '-'■*' /o-jj^ = 
Y(»)/[(<j«)2 _^ y(^))2] ^j^j^ r^{,) ^ ^(«)/(^W)2^ leading to 



where we introduced the current as j^'^' = a. — — which 
is constant through the whole chain. Although this equa- 
tion is equivalent to Eq. (|l|) it has a different physical in- 
terpretation, since it describes the heat flux through the 
granular array. In the following we consider the case of 
constant current. 



A. Kinetic heat equation in dimensionless units 

It is convenient to rewrite the kinetic Eq. (Ef) in dimen- 
sionless units. Temperatures and voltages are measured 
in units of the first lead temperature Ti and the time 
scale is in units of mean energy level spacing (5, i.e. 

yw = ey«/(fcBri), T = ts/h, (5) 

where Tg , T2, "Is are the grain, second lead, and sub- 
strate dimensionless temperatures, t the dimensionless 
time, and V*-*-' the dimensionless voltage drops across 
grain i. We also write the conductivities, Eqs. (g), as 
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where the 'tilde' expressions are dimensionless. The di- 
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and wc define the dimensionless thermal conductivity as 
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where we introduced the dimensionless parameters' 
2 / m \ 2 






/3 = 37ph 



a 



ph 



Substituting Eqs. (|6|) back into Eq. ( 
following dimensionless heat equation 
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where V is the dimensionless gradient and the dimen- 
sionless electron-phonon coupling term to the substrate is 
5 
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Equation (Bh is the main equation which is considered in 
the following. 
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The total voltage drop over the chain is then Vtot = 
Y^^ V^'^. The l.h.s of Eq. (||) is explicitly given by (writ- 
ten as central differences; cf. also Eq. (23) of Ref. [||) 
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Using Eqs. (^0|) and (0) we obtain for the gradient of 
thermal conductivity Vk'*' the following expression 
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Here we took into account the fact that the tunneling 
conductance g^ is an averaged quantity for a particular 
grain (i), with 
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where V*^*^ is given by Eq. (M. Using Eqs. (|l])-(|3|) we 
solve the kinetic heat Eq. ^ numerically for different 
situations. In the simulations we keep the lead temper- 
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atures Ti = Tg 



1 and T2 = Tl 
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fixed and inte- 



grate (pf) until a steady state is reached. For the initial 
condition of the grain temperatures we choose a linear 
interpolation from Ti to T2. 



C. Kinetic heat equation: Stationary case 

Before solving kinetic heat Eq. (m for the dimension- 
less temperatures Tg (r) vs. time r exactly, we dis- 
cuss its stationary solution corresponding to large times, 
T ^ T* . where r* is some characteristic time scale which 



we will discuss in section [II E , We also consider only the 
case *Ps ~ (no coupling between the array and the 
substrate) in this section. In other words we first esti- 
mate how pronounced the effect of inelastic cotunneling 



on grain heating is. For large times the grain tempera- 
ture reaches some constant steady state value, 



lim s:«(t) = s:«* 
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i.e. d'Tg /dr = 0. In this section we consider only 
the stationary case and suppress the superscript * in the 
following. Using Eq. (H) we obtain 
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In the limiting cases of low voltage (V^*-* ^ 'Jg ) or low 
temperature (iL ^ V'*-*) one can find the analytical so- 



lution of Eq. 

r(') 



for the dimensionless grain temperature 
Tg'''. However, due to the cumbersome expression for the 



discrete heat-flux, Eq. (nil) we go over to the straight 
forward continuum approximation, i.e. the grain index 
i goes over in coordinate x and the discrete gradients 
become continuum ones. 



1. Low voltages 

We first consider the stationary solution, Eq. 
the limit of small voltages, V <C Tg. Using Eqs. 
and (p^ ) we obtain 
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FIG. 3. (Color online) Stationary heat profiles (grain temper- 
ature Tg vs. grain index) for Ti = Ts = lOK for asymmetric 
case with different right lead temperature T2, in particular 
T2 > Ti. The curves for tunneling conductance gt — 0.5 for 
different chain lengths without substrate coupling (*Ps — 0) 
are chosen in such a way that the right-lead temperature T2 
lies on the symmetric profile with T2 — lOK (number of grains 
iV = 100), i.e. T2 = 12.4K for iV = 50 and T2 = 11. 8K for 
TV = 80. For the latter two curves the end is marked by a 
vertical line since all three curve lie on top of each other. The 
second set of three curves is plotted for gt =0.1 and electron- 
phonon coupling constant Es — 10*W/(ni"'K^), for three dif- 
ferent right lead temperatures T2 = 10, 12, 14K, written next 
to these curves (all for A^ — 100). The current density is fixed 
at j — lA/cm^ for all graphs. 
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FIG. 2. (Color online) Stationary heat profiles (grain tem- 
perature Tg vs. grain index) for the symmetric case where 
both lead temperatures coincide Ti — T2 — lOK and for dif- 
ferent tunneling conductances gt — 0.5,0.1,0.05. The graph 
for gt =0.1 is almost identical with the profile for random 
gt with mean 0.1 and variance 0.05 (dotted, magenta) and 
the approximation for gt = 0.1 (dash-dot, cyan). A curve 
with enabled coupling to the substrate (see definition below 
Eq. (|)) is shown as well for gt = 0.1 and E^ = 10''W/(m^K^) 
at a substrate temperature of Ts = lOK. The current density 
is fixed at j — lA/cm^ and length A'^ = 100. 



where we introduce the notation T? = .>, ■> ^'^ , „ ; ^ . 



Equation (16) can be integrated using the new variable 
z(Tg) ^ I^VTg. Then Vz = z'VTo 



■g ~ z^z/T^, where 
z' — dz/df£g. In terms of variable z, Eq. (|l6| ) has the 
form 
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leading to the solution z = T/,, /T^„ — T^, with T™ being 

the integration constant, physically meaning max^, Tg(.T) 
which has the coordinate Xm, i.e. Tm = Tg(a;m). Using 
this solution we obtain the following differential equation 
for the dimensionless temperature 
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leads to the following solution 

1 
3 

Introducing the new variable y — Tg/T^, Eq. ( |l9| ) can 
be written as follows 

y3 + 3y2 = 4 - (x - x™)2(3Ti/T3j2 = s + 2. (20) 



The positive root of Eq. ( pO| ) reads 
2/3 = ^1/3 (s) - 1 



(21) 



where we introduced Chebyshev's cubic root function, 
see Appendix ^ 
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Using Eq. (^) for the dimensionless grain temperature 
we finally obtain 



'Igix) 



^JCi/3 (s) - 1, 



(23) 



where parameter T^ is defined below Eq. (Kw. For the 
symmetric case (both lead temperatures are equal, Ti — 
T2) Eq. (p3|) is plotted in Fig. 0. One can see that the 
maximum temperature Tm is reached in the middle of 
the chain. 



2. High voltages 

Here we discuss the stationary solution, Eq. (iq), in 
the limit of small temperatures, Tg ^ V. Using Eqs. (0) 
and (nS) we obtain the following equation 
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where we introduce the notation T'Jr = — tttt. 

tion (El) can be integrated using the new variable. 



Z{1g) = IgVlg. 

z' = dz/dfZg. 
form 



Then Vz = z'Vlg = z'zj^g, where 
In terms of variable z, Eq. (p4) has the 
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leading to the solution z = T^^/'X^ — 1^. Equation ( |25| ) 

is similar to Eq. ( |l7| ) obtained in the large temperature 
limit. The solution of Eq. (|25|) reads 



1g{x) 



with s 



(x-x„02(3Tff/T3j 



Cl/3 (5) - 1. 



3 \2 



(26) 



> 

o 










N=100 
5!05 

=10%/ 


Ti=Ts= 


10K, 


j=1A/cm^ 


9r 


=0.1 


9t= 
2s 


- 








91=0-1" 
g,=o.5 


— 


:=:: 



0.25 
0.2 

> 0.15 

O 

>" 0.1 

0.05 


10 15 20 25 30 35 40 45 50 

T2[K] 

FIG. 4. (Color online) Voltage dependence of nano dot chain 
Vtot vs substrate temperature Ts (top) and right lead tem- 
perature T2 (bottom) for different tunneling conductances 
Qt- All plots are for applied current density j = lA/cm^, 
left lead temperature Ti = lOK, and chain length A*' = 100. 
In the upper panel the electron-phonon coupling constant is 
Es = 10^W/(m^K^). In the lower plot the electron-phonon 
coupling is switched off, apart from the dotted blue curve 
marked with gt = 0.1, E^ = 10^W/(m^K^) which is coupled 
to a substrate at Ts = lOK. 



is slightly above the lead temperature. It is clear that the 
sensitivity is reduced if only one lead is used as tempera- 
ture sensor. In the case of the substrate sensor, the sen- 
sitivity range can be tuned by the tunneling conductance 
of the chain. Another important issue for applications as 
temperature sensors is the response time of the device. 
In the case of the chain of grains, this response time can 
be reduced by (i) increasing the current density, j, (ii) 
having a better electron-phonon coupling to the sample, 
i.e. larger Sg, (iii) or smaller tunneling conductances, gt] 
see next section. Timescales of several nanoseconds are 
typical for not too long chains. The length of the chain 
can be adjusted for best sensitivity and fastest response 
time. 



D. Temperature-voltage dependence 



As mentioned in the introduction the chains of nano 
dots can be used as highly sensitive thermometers from 
cryogenic up to room temperature. In order to demon- 
strate this, we studied the total voltage drop, Vtot, on the 
chain for different substrate Tg or lead T2 temperatures. 
The numerical results for Vtot are shown in Fig. H for Ts 
and T2 in the range from lOK to 50K. As one can see, the 
chains are particularly sensitive if the substrate (sample) 



E. Kinetic heat equation: Dynamic case 

In this section we estimate the characteristic dimen- 
sionless timescale t* [see Eq. (@)] associated with the 
heating of nano dot chain due to inelastic cotunneling. 
We define t* as the typical timescale to reach the max- 
imum temperature T^ in the chain, 1^ = max^ Tg = 
linv^oo ^^m (''■)• The time r* can be easily estimated for 
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FIG. 5. (Color online) Plot of the rescaled maximum grain 
temperature Tm{t) vs. time t. The typical timescale t* to 
reach the steady state value T^ is shown for tunneling con- 
ductance gt — 0.5 and current density j = lA/cm^ (solid, 
red curve). For comparison also the graphs for two different 
substrate coupling are shown: for larger coupling the typical 
timescale t* is reduced. The plot also shows the time evolution 
of Tm{t) for different gt and j (see legend) for zero electron- 
phonon coupling (no substrate): for smaller gt or larger j, the 
timescale f* can be decreased as well. 
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FIG. 6. (Color online) Upper panel: Solution of Eq. (G3) for 
maximum temperature T^ vs. lengths N of nano dot chain 
for a fixed tunneling conductance gt = 0.1. Low panel: T^ 
vs gt for a fixed length A'' = 100. Both graphs are plotted for 
constant current density j — lA/cm^. 



low voltages V <C T. Using the fact that for symmetric 
case (both leads have the same temperature, Ti = ^2) 
the maximum temperature Tm is located in the middle 
of the chain {i = N/2), meaning that V'Xm = 0, and 
using Eq. (|^) in the absence of a substrate, *Ps = 0, we 
obtain the following equation 
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Integrating Eq. ( p7| ) for a given fixed current t and us- 
ing Eq. (|6|) for the dimensionless conductivity CTj^^ we 
obtain the following result 
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Here the maximum temperature T*j is implicitly defined 
by the following equation [see also Eq. (P3|)] 
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The solution of Eq. (^, TJ^ vs number of grains, N, and 
TJ^ vs. tunneling conductance, gt, is shown in Fig. ^ 



IV. OVER-HEATING HYSTERESIS 

If the rate of the energy supply from the external bias 
to the charge carriers exceeds the rate of energy losses to 



the environment, the phenomenon of overheating appears 
and the energy distribution function of the current carri- 
ers noticeably deviates from the equilibrium distribution 
functionJ23'E3 One of the characteristic manifestations of 
the overheating effect is the onset of the "falling" re- 
gion of the I-V curve where the differential conductivity 
G = dl/dV < 0. The corresponding I{V) character- 
istics is referred to as that of the S'-type if the current 
is the multi-valued function of the voltage, and the I- 
V curve of the A'^-type corresponding to the case where 
the current is nonmonotonic but still remains a single 
valued function of the voltage. The phenomenon of over- 
heating has been a subject of the incremental interest 
and extensive studi^ during the decades, see review of 
Volkov and KoganJ23 and the impressive progress in un- 
derstanding of the underlying mechanisms was achieved. 
Recently it was suggestedcj that the mechanism of heat- 
ing instability described in detail in Ref. E3 combined 
with the expression for the dissinatci power calculated 
for strongly disordered conductorsEErO (see Ref. g for re- 
view) can explain switching I-V characteristics in InO 
and TiN samples. 

Using Eq. (g|) we simulate the hystcrcstic overheating 
behavior in chains of nano-grain arrays. We first increase 
the current from j = OA/cm^ to 5A/cm^ and then de- 
crease it back to zero. The results are shown in Fig. M. 
One can see a pronounced hysteretic behavior of both the 
total voltage drop Vtot and the maximum chain temper- 
ature T^. It is important to remark that this hysteretic 
behavior disappears for shorter chains. 



N=100 




[A/cm^] 



FIG. 7. (Color online) Hysteretic behavior of the maximum 
temperature r,*^ (upper panel) and the total voltage trop Kot 
(lower panel) vs current density j. The current density j was 
first increased from OA/cm^ to 5A/cm^ and then decreased 
back to OA/cm^ . The arrows show the direction of the hys- 
teresis. The red curves correspond to the case with no cou- 
pling to the substrate. The green curves are shown in the 
presence of the substrate, with electron-phon coupling con- 
stant Es = 10^W/(m^K^). The tunneling conductance for all 
curves was fixed, gt — 0.1. 



V. DISCUSSIONS 

In all simulations throughout this work, we use typical 
experimental parameters to estimate the heating effects 
on a chain of weakly coupled grains and presented all 
equations in either SI or dimensionless units. In par- 
ticular we used the following physical parameters: the 
grain size a = lOnm, Debye temperature Od ~ 450K, 
Coulomb energy Ec « 1600K (with relative static per- 
mittivity K = 1), mean energy level spacing S = IK, 
phonon mean free path Zph = a, and tunneling conduc- 
tance gt € [0.05; 0.5]. Furthermore we simulated typical 
chains of A^ = 100 grains (i.e. of a few fira total length). 

For the electron-phonon coupling to the substrate we 
used values of 10^ or 10^W/(m'^K^) since we want to de- 
scribe only the weak coupling of the chain to a different 
substrate. For comparison, typical values for bulk ma- 
terials are much larger (on the order 10^W/(m'^K^), see 
e.g. i). 

As for the analytical approximations in the low- and 
high-voltage regions (without substrate coupling) the 
particular form of the differential and later algebraic 
equations, allowed us to derive a closed solution using 



Chcbyshev's cube root. 

As a first step we calculated the stationary heat profile 
of the chain under different conditions. It is important to 
remark that the low-voltage approximation fits very well 
the full simulation and randomness in the tunneling con- 
ductance has only minor effects on the final profile (pro- 
vided the chain is not "broken" by a very small tunneling 
conductance.) Furthermore we estimated the character- 
istic heating timescale t* and maximum temperature T,^ 
in Eq. (Eq) using the analytical solution and obtained as 
a result t* « 0.1/xs and T^^ « 20K. 

These estimates can be directly compared with our 
simulations of Eq. (Bf) for rescaled maximum grain tem- 
perature Tm{t) VS. time t. The outcome of these sim- 
ulations is shown in Fig. H. One can see that the typ- 
ical timescale t* to reach the steady state value T*-^ is 
of order of 0.1/xs. For completeness in Fig. || we also 
presented several different curves Tm{t) vs. time t for 
different tunneling conductances gt and current densi- 
ties j. In addition, we also presented two curves in the 
presence of the substrate: for larger coupling the typi- 
cal timescale t* is reduced. Based on the above there 
is a good agreement between our estimates and simula- 
tions. At this point we note that in general the kinetic 
heat Eq. (g) has an additional term proportional to the 
thermoelectric coefficient. Recently, thermoelectric Ppfi>Pr 
erties of weakly coupled granular arrays were studiedEHlZI. 
In particular, in Ref. ]13 we estimated the thermoelectric 
and Seebeck coefficients. It was found that the thermo- 
electric coefficient has an additional small factor propor- 
tional to T/Ep <C 1 compared to the electric and thermal 
conductivities, where Ep is the Fermi energy, i.e., for the 
temperatures under consideration this additional term in 
Eq. (|9|) would be about three orders of magnitude smaller 
than the others. 

In conclusion, we have studied heating effects in one- 
dimensional chains of quantum dots due to inelastic elec- 
tron cotunneling using a kinetic heat equation. In the 
low and high voltage limits we solved the stationary heat 
equation analytically. We demonstrated the possible ap- 
plication of chains of nanograins as highly sensitive ther- 
mometers and estimated the typical timescale to reach a 
steady heat profile in these chains. Finally, we showed 
the over-heating hysteresis in the large-current or voltage 
regimes. The infiuence of the electron-phonon coupling of 
the chain to a substrate was considered numerically and 
has pronounced effects on the stationary heat profile, the 
typical timescales to reach it, and the thermometric and 
hysteretic behavior. 
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some. In this paper [see Eq. (^0|)] we can remark ably 
find a simple expression for the root because Eq. ( Al) 
reduces to the problem of finding roots of the Chebyshev 
polynomial. T he tr ansformation y = y — 1 reduces the 
left side of Eq. ( Al) to the third Chebyshev polynomial, 



Cs^y) ^ y — 3y, and the right side to s. Then the re- 
sulting equation has only one positive root known as the 
Chebyshev cube root:E£l y = Ci/3(s), where 



Ci/3(s) = 2 cosh 



- acosh(s/2) 
o 



(A2) 



FIG. 8. (Color online) The graph of Ci/3(s) (solid line). The 
graphs with the dashed and dotted lines sh ow t he behavior of 
the asymptotical approximations, see Eq. (A3). 



Appendix A: Cube roots 

Here we discuss the positive roots of the cubic equation 



y^ + 3y'-^s + 2. 



(Al) 



The roots of the cubic equation can be always found 
analytically however the expressions are quite cumber- 



Here acosh(s/2) ~ ln[(.s^ -I- Vs^ — 4)/2] and we use the 
branch of the logarithm which is real on the positive real 
line and the branch of the square root which is positive 
on the real axis. Asymptotic behavior of the Chebyshev 
cube root is the following: 



^1/3 (s) 



1+V^' if0<s + 2<l; 
2cosh(i|^), s> 1. 



(A3) 



The graph of Cit^ is shown in Fig. 
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